GLOBAL WELL-POSEDNESS FOR EULER-BOUSSINESQ 
SYSTEM WITH CRITICAL DISSIPATION 



T. HMIDI, S. KERAANI, AND F. ROUSSET 

Abstract. In this paper we study a fractional diffusion Boussinesq model which 
couples the incompressible Euler equation for the velocity and a transport equa- 
tion with fractional diffusion for the temperature. We prove global well-posedness 
results. 



1. Introduction 

Boussinesq systems of the type 

dtv + v • Vv + Vp = 9e2 + vT> v v 
dtO + v-Ve = kVqQ 
div v = 

v lt=0 = v°, e lt=0 = e° 

are simple models widely used in the modelling of oceanic and atmospheric motions. 
These models also appear in many other physical problems, we refer for instance 
to [3 [3] for more details. Here, we focus on the two-dimensional case, the space 
variable x = (xi,X2) is in M 2 , the velocity field v is given by v = (v , v ) and the 
pressure p and the temperature 9 are scalar functions. The factor 9e2 in the velocity 
equation, the vector e<i being given by (0, 1), models the effect of gravity on the 
fluid motion. The operator T> v and T>g whose form may vary are used to take into 
account the possible effects of diffusion and dissipation in the fluid motion, thus the 
constants v > 0, k > can be seen as the inverse of Reynolds numbers. 
Mathematically, the simplest model to study is the fully viscous model when v > 0, 
K > and T> v = A, Dg = A. The properties of the system are very similar to the 
one of the two-dimensional Navier-Stokes equation and similar global well-posedness 
results can be obtained. 

The most difficult model for the mathematical study is the inviscid one, i.e. when 
v = k = 0. A local existence result of smooth solution can be proven as for symmetric 
hyperbolic quasilinear systems, nevertheless, it is not known if smooth solutions can 
develop singularities in finite time. Indeed, the temperature 9 is the solution of a 
transport equation and the vorticity uj = curl v = d\v 2 — &2V 1 solves the equation 

(1.1) d t u + v ■ = d x 9. 
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The main difficulty is that to get an L°° estimate on u which is crucial to prove 
global existence of smooth solutions for Euler type equation, one needs to estimate 
Jo ll^i^ll-L 00 an d, unfortunately, no a priori estimate on d\9 is known. 
In order to understand the coupling between the two equations in Boussinesq type 
systems, there have been many recent works studying Boussinesq systems with par- 
tial viscosity i.e. with a viscous term only in one equation. For k > 0, v = and 
T>0 = A, the question of global existence is solved recently in a series of papers. In [7j, 
Chae proved the global existence and uniqueness for initial data (t>°, 0°) G H s x H s , 
with s > 2, see also |20| . This result was recently extended in |16| by the two first 

-+i — 1+- 

authors to initial data v G B pl and 9° G B pl p n L r ,r > 2. More recently, the 
study of global existence of Yudovich solutions for this system has been done in |14j . 
We also mention that in [15] , Danchin and Paicu were able to construct global strong 
solutions (still for k > 0, v = 0) for a dissipative term of the form T>q = d\\9 instead 
of A9. Recently the first author and Zerguine [19] proved the global well-posedness 
for fractional diffusion T>q = — |D| a for a G]l,2[ where the operator |D| a is defined 
by 

H\Wu){i) = \Z\ a (Fu)(0. 

In these works, the global existence result relies on the fact that the only smoothing 
effect due to the transport- (fractional) diffusion equation 

d t + v-V0 + \B\ a = O 

governing the temperature is sufficient to counterbalance the amplification of the 
vorticity. However the case a = 1 is not reached by their method. The main reason 
is that this case can be seen as critical in the previous approaches in the sense that 
the smoothing effect for the temperature equation does not provide the Lj,(L°°) 
bound for d\9 which seems needed to control the amplification of the L°° norm of 
the vorticity. Note that such an estimate is nevertheless almost true since d\9 can 
be estimated in the space L^(L°°) which has the same scaling (see below for the 
definition). The aim of this paper is the study of the well-posedness for this case, 
i.e., we focus on the system 

dtv + v ■ Vv + Vp = 9e2 
d t 9 + v-V9+\B\9 = 
divv = 

v ]t=0 = v°, 9 {t=0 = 9°. 

Note that we have taken k = 1 which is legitimate since we study global well- 
posedness issues. Indeed, we can always change the coefficient n into 1 by a change 
of scale. 

We also point out that at first sight, the system (jl.2p contains the mathematical 
difficulties of the critical quasigeostrophic equation introduced in [10] 



(1.2) 



(1.3) d t 9 + v V0 + |D|0 = O, v = V ± \D\~ 1 9 

which was much studied recently. Indeed, in (|1.2j) the link between v and 9 is not 
given by the Riesz transform but by a dynamical equation, the first equation of 
(|1.2j) . Nevertheless, from this velocity equation one gets that v has basically the 
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regularity of 9 as in the quasigeostrophic equation. The global well-posedness for 
(|1.3p was obtained recently by Kiselev, Nasarov and Volberg |21j . We also refer to 
the work [6] by Caffarelli and Vasseur about the regularity of weak solutions. Other 
discussions can be found in [T| fill [12] . 

The main result of this paper is a global well-posedness result for the system (|1.2p 
(see section [2] for the definitions and the basic properties of Besov spaces) . 

Theorem 1.1. Let p g]2,oo[, v° G Bl ol PiW 1,p be a divergence- free vector field 
ofM? and 9° G B^ 1 n L p . Then there exists a unique global solution (v,9) to the 
system (jl.2p with 

v G L£ c (R+; Bl A n W x *\ 9 G L£ C (R+; B^ n LP) n ^(M+j B^). 
A few remarks are in order. 

Remark 1.2. If we take 9 = then the system (jl.2p is reduced to the well-known 
2D incompressible Euler system. It is well known that this system is globally well- 
posed in H s for s > 2. The main argument for globalization is the BKM criterion 
[2] ensuring that the development of finite time singularities for Kato's solutions 
is related to the blowup of the L°° norm of the vorticity near the maximal time 
existence. In [26j Vishik has extended the global existence of strong solutions result 

to initial data lying on the spaces 

B l^i /P - Notice that these spaces have the same 
scaling as Lipschitz functions (the space which is relevant for the hyperbolic theory) 
and in this sense they are called critical. We emphasize that the application of the 
BKM criterion requires a super- lipschitzian regularity (H s with s > 2 for example). 

For that reason the question of global existence in the critical spaces B^^ V is hard 
to deal with because these spaces have only a lipschitzian regularity and the BKM 
criterion cannot be used. 

Since Bt + ^ p B^^ n W 1,r for all p G [1, +oo[ and r > max{j>, 2}, then the space 
of initial velocity in our theorem contains all the critical spaces B^ 2 ^ except the 
biggest one, that is B^ 1 . For the limiting case we have been able to prove the global 
existence only up to the extra assumption Vu° G LP for some p G]2, oo[. The reason 
behind this extra assumption is the fact that to obtain a global L°° bound for the 
vorticity we need before to establish an LP estimate for some p G]2, oo[ and it is not 
clear how to get rid of this condition. 

Remark 1.3. Since Vt>,V# G L 1 1 oc (IR_|_; L°°) (see Remark 15.61 below for 9) then we 
can easily propagate all the higher regularities: critical (i.e. vq G B^ 2 ^ with p 
finite) and sub-critical (for example vq G H s , for s > 2). 

The main idea in the proof of Theorem ll.il is to really used the structural properties 
of the system solved by (w, 9), lo = curl v = d\v 2 — 82V 1 . Indeed, if we neglect the 
nonlinear terms for the moment, one gets the system 

8^ = 8x9, 8 t 9 = -\D\9 
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and we notice that its symbol given by 




is diagonalizable for £ ^ with two real distinct eigenvalues which are and — 1£|. 
By using the Riesz transform 1Z = <9i/|D|, one gets that the diagonal form of the 
system is given by 

d t ne = \B\ne, dt{u + ne) =0. 

This last form of the system is much more convenient in order to perform a priori 
estimates. To prove Theorem ll.il we shall use the same idea, we shall diagonalize the 
linear part of the system and then get a priori estimates from the study of the new 
system. The main technical difficulty in this program when one takes the nonlinear 
terms into account is to evaluate in a sufficiently sharp way the commutator \TZ, v ■ 
V] between the Riesz transform and the convection operator. Such commutator 
estimates are stated and proven in section [3] of the paper. 

The diagonalization approach used in this paper also allows to prove global well- 
posedness in different spaces for a "Boussinesq-Navier-Stokes" system i.e. the system 
which corresponds to Vg = and V v = — |D|. This is discussed in a companion 
paper [18] . 

The remaining of the paper is organized as follows. In section [2] we recall some 
functional spaces and we give some of their useful properties. Section [3] is devoted 
to the study of some commutators involving the Riesz transform. In section |4] we 
study a linear transport- (fractional) diffusion equation. Especially, we establish 
some smoothing effects and a logarithmic estimate type. 

In section [5] we give the proof of Theorem 11.11 which is It is splitted into three parts. 
We first establish some a priori estimates, then we prove the uniqueness and finally 
we briefly explain how one can easily combine a procedure of smoothing out of the 
initial data with the a priori estimates to get the existence part of the theorem. An 
appendix is devoted to the proof of a technical commutator lemma. 

2. Notations and preliminaries 

2.1. Notations. Throughout this work we will use the following notations. 

• For any positive A and B the notation A < B means that there exist a positive 
harmless constant C such that A < CB. 

• For any tempered distribution u both u and J-u denote the Fourier transform of 
u. 

• Pour every p £ [l,oo], || • \\u> denotes the norm in the Lebesgue space L p . 

• The norm in the mixed space time Lebesgue space LP([0, T],L r (R d ) is denoted by 
II ' Hz4L r (with the obvious generalization to || • \\u^x f° r an y normed space X). 

• For any pair of operators P and Q on some Banach space X, the commutator 
[P, Q] is given by PQ - QP. 

• For p £ [1, oo], we denote by W 1,p the space of distributions u such that Vu € BP . 
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2.2. Functional spaces. Let us introduce the so-called Littlewood-Paley decom- 
position and the corresponding cut-off operators. There exists two radial positive 
functions % 6 and 93 G £>(R d \{0}) such that 

i) X(0 + Yl ¥>(2" 9 = !; V g > 1, supp X n supp ^(2" 9 ) = 

<?>o 

ii) supp (p{2~j-) n supp ^(2~ fc -) = 0, if |j - fe| > 2. 
For every u G ^'(M^) we set 

A_it; = x(D)f ; Vg G N, A g w = ip(2~ q B)v and 5, = ^ Aj. 

i=-i 

The homogeneous operators are defined by 

A q v = ip{2~ q B)v, S q v = ^ Ajv, \/q G Z. 

From [5] we split the product uv into three parts: 

uv = T u v + T„it + R(u, v), 

with 

l 

T u v = ^Sq-iuAgV, R(u,v) = ^ A q uK q v and A g = Aq +i . 
g g i=-l 

For (p, r) G [l,+oo] 2 and sGKwe define the inhomogeneous Besov space 5® as 
the set of tempered distributions u such that 

lldls* := (2 qs \\A„u\\ L p ] < +oo. 



The homogeneous Besov space fi® is defined as the set of u G 5'(R ) up to poly- 
nomials such that 



U c 



2 9S ||A g M|| L p < +oo. 



P,r \ / Br 

Let T > and p > 1, we denote by Lj.B^ r the space of distributions it such that 



^ s ||A 3 n||iP 
We say that it belongs to the space L^IBp r if 



T p,r 

< +0O. 



IMIZ^ r : = { 2qS W^A\LP T Lv) r < +00. 

By a direct application of the Minkowski inequality, we have the following links 
between these spaces. Let e > 0, then 

L^,Bp r c — ► L^,Bp r c — > L^Bp r £ , if r > p, 

Lq^Bp^ 6 c — ► LqiBp r c — ► Lq^Bp r , if p v. 
We will make continuous use of Bernstein inequalities (see [8] for instance). 
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Lemma 2.1. There exists a constant C such that for q,k £ N, 1 < a < b and for 

f G L a (R d ), 

sup \\d a S q f\\ Lb < C k 2^ k+d ^-D)\\S q f\\ L a, 

\a\=k 

C~ k 2 qk \\/\ q f\\ La < sup \\d a A g f\\ L a < C k 2 qk \\k q f\\ La . 

\ct\=k 

3. RlESZ TRANSFORM AND COMMUTATORS 

In the next proposition we gather some properties of the Riez operator 1Z = d\/\D\. 

Proposition 3.1. Let 1Z be the Riez operator 1Z = d\/\D\. Then the following hold 
true. 

(1) For every p e]l, +oo[, 

11^ II £(lp) ~ !■ 

(2) Let x G T>(R d ). Then, there exists C > such that 

|||D| S X (2- 9 |D|)^|| £(LP) <C2^, 

for every (p,s,q) € [1, oo] x]0, +oo[xN. 

We notice that the previous results hold true if we change |D| S by V s 
with seN. 

(3) Let C be a fixed ring. Then, there exists ip £ S whose spectum does not meet 
the origin such that 

TZf = 2 qd i){2 q -)*f 
for every f with Fourier transform supported in 2 q C. 

Proof. (1) It is a classical Calderon-Zygmund theorem (see |24j for instance). 
(2) Let K G S' such that FK{£) = |£| s x(l£l)fj- K is a tempered distribution such 
that its Fourier transform TK is C 00 (M <:i \ {0}) and satisfies for every a £ N d 

\d^K(0\<C a \Z\ s -^. 

According to Mikhlin-Hormander Theorem (see [21] for instance) we have for every 
a G N d and x G R d \{0}, 

\d%K(x)\ < C' a \x\- S ~ d ~^. 

Then, it ensues that 

\K(x)\ < C\x\- d ~ s 

Since J-K is obviously in L 1 , we also have that Li G Co(M ). This removes the 
singularity at the origin and gives 

\K(x)\ <C(l + \x\)- d - s Vx G M. d . 

Therefore we get that the kernel K G L . Now for every u G TP we have |D| s 7?-x(|D|)n 
K -ku. We can now conclude the case q = by using the classical Young inequality 
for convolution products. 

The case q > 1 can be derived from q = via an obvious argument of homogeneity. 
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(3) This is can be done easily by introducing a judiciously chosen cut-off function. □ 



The following lemma will be useful in the proof of many commutator estimates. 

Lemma 3.2. Given (p,m) G [l,oo] 2 such that p > m! with m' the conjugate ex- 
ponent of m. Let f,g and h be three functions such that V/ G L p ,g G L rn and 
xh G L m ' . Then, 

\\h*(fg) - f{h*g)\\LP < \\xh\\ Lm > \\Vf\\ LP \\g\\ L m. 
Proof. We have by definition and Taylor formula 

h*(fg)(x)-f(h*g)(x) = [ h(x-y)g(y){f(y)-f(x))dy 

= g{y)K x ~ v) (V ~ x ) -V/(a; + t(y - x)) dydt. 

Jo JR d L J 

Using Holder inequalities and making a change of variables z = t(x — y) we get 

"(/ 

where we set h\{z) = \z\\h(z)\. Using Convolution inequalities we obtain since p > m' 



h*(fg)(x)-f(h*g)(x)\ < \\g\\ Lm f ( / t' d hf ' {t~ l z)\V f\ m ' (x - z)dz) 



\h*(fg)-f(h*g)\\ LP < Nl^ll^ll^lllV/ni^ 

L m' 

< \\9\\l*" \\hl || rm' || V/||l>. 



L" 

□ 

As explained in the introduction, the control of the commutator between 1Z and the 
convection operator v ■ V is a crucial ingredient in the proof of Theorem 11.11 

Theorem 3.3. Let v be is a smooth divergence- free vector field. 

(1) For every (p,r) G [2, oo[x [1, oo] there exists a constant C = C(p,r) such that 

m,v ■ vn B o r < cnvt>Mii0|| Bgoir + ii^iilp), 

for every smooth scalar function 9. 

(2) For every (r, p) G [1, oo] x]l, oo[ and e > there exists a constant C = 
C(r,p,e) such that 

\m,v v]9\\ B o xr < ccimIloo + \\u\\ LP )(\\e\\ B ^ r + \\e\\ LP ), 

for every smooth scalar function 9. 
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Proof. We split the commutator into three parts, according to Bony's decomposition, 

[n,v-v]e = ^[K.Vuj.vjV + ^^V.viVifl 

qGN q&N 

+ ^[K,A q vV]A q e 

q>-\ 

= E^ + E 11 ^ E m « 

gGN (jSN q>-l 

= I + 11 + III. 

We start with the estimate of the first term I. According to the point (3) of 
Proposition 13.11 there exists h £ S whose spectum does not meet the origin such 
that 

l q (x) = h q * (S q -!V ■ VAqO) - Sg-xv ■ (hg * VA q 0), 
where h q (x) = 2 dq h(2 q x). Applying Lemma 13.21 with m = oo we get 

(3-1) < ||Vt»|| LP ||A^|| LO c. 

In the last line we've used Bernstein inequality and ||x/t ? ||ii = 2 _9 ||x/i||£i . 
Combined with the trivial fact 

^Ei,= E l * 

i b'-?|<4 

this yields 

L q\\LP 

q>-l 

^ II II^IIb^. 

Let us move to the second term II. As before one writes 

llg(x) = hg * {AgV • VSg-tf) ~ AgV • (hg * V S ' g-rf) , 

and then we obtain the estimate 

||ILJ LP < 2~' } \\AgVv\\ L p\\Sg- 1 Ve\\ Laa 
< \\Vv\\ LP 2^^||A^|| L o . 

Combined with convolution inequalities this yields 

1 1 II 1 1 < ||Vu |Up||0||_bo . 

11 p,r oo,r 

Let us now deal with the third term III. Using that the divergence of A q v vanishes, 
we rewrite III as 

III = ^^div(A g uA^)-^div(A 9 u^A^) + ^[^,A g u- V]A g 

q>2 q>2 q<l 

= Jl + J 2 + J 3 - 



\ i \\b^ < ( E iiv 
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Using Proposition 13. 11 - (2) . we get 

\\Aj1Zdiv(A q v A q 9)\\ LP < 2 J '||A^|| LP ||Ag0|| L o O . 
Also, since A q 9 is supported away from zero for q > 2 then Proposition [37T] (3) yields 

|| Ajdiv(A 9 t> TZA q 9)\\ LP < 2 j \\A q v\\ LP \\KA q 6\\ L oo 

< 2 j \\A q v\\ LP \\A q e\\ L ~. 

Therefore we get 

\\A j (J 1 + J 2 )\\ L p < Y 2 j '||A 9 d|| Lp ||A^|| Z/ ^ 

q£N 
V>3— 4 

< \\Vv\\ L v Y 2 J - 9 ||A^|| LO c, 

q£N 
q>j-4 

where we have again used Bernstein inequality to get the last line. It suffices now 
to use convolution inequalities to get 

||«7i + ^2||bo < I|Vi;||lj>||0||bo . 

For the last term J3 we can write 

Y [K, \ v ■ V]A^(x) = ^[div x(D)^, A q v]A q e(x), 
-i<?<i g<i 

where x belongs to T>(M, d ). Proposition 13. II ensures that div x{T))lZ is a convolution 
operator with a kernel h satisfying 

\h(x)\<(i+\x\r d -\ 

Thus 

J 3 = Y h * (\ V " - A q V ■ (Jl * A q 0). 

First of all we point out that A, J3 = for j > 6, thus we just need to estimate the 
low frequencies of J3. Noticing that xh belongs to LP' for p' > 1 then using Lemma 
13.21 with m = p > 2 we obtain 

I|AjJ 3 ||l°° < Y \\ x M\lp' \\A q Vv\\ L p\\A q e\\ LP 

^ ||Vu||£p ||A,0||i*. 
-i<?<i 

This yields finally 

ll^llso^ < ||Vu||lp||0||lp. 

This completes the proof of the first part of Theorem 13.31 The second part can be 
done in the same way so that we will only give here a shorten proof. To estimate 
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the terms I and II we use two facts: the first one is ||AgVit||i,°o ~ ||A g a;||ioo for all 
q £ N. The second one is 

q-2 

WVS g -iv\\L^ < ||VA_it>||x°° + ^ HAjVvHi.oo 

< \\u)\\lp + g||w||L°°- 

For the remainder term we do strictly the same analysis as before except for J3: we 
apply Lemma 13.21 with p = 00 and m = p leading to 

WAjMIlp < y~) \\xh\\ LP > 11 A g V«||L°° i| A q e\\ LP 

9<1 

^ ||Vw||iP \\ A <i e h p 

-i<g<i 
< ||w||lp||^||lp. 



This ends the proof of the theorem. 



□ 



4. Transport-Diffusion equation 



In this section we will give some useful estimates for any smooth solution of a linear 
transport-diffusion model given by 



(TD) 



d t 8 + v -X79 + \D\8 = f 
0, t=o = 9°. 



We will discuss three kinds of estimates: LP estimates, smoothing effects and loga- 
rithmic estimates. 

The proof of the following LP estimates can be found in |13j . 

Proposition 4.1. Let v be a smooth divergence-free vector field ofM. d and 6 be a 
smooth solution of (TD). Then we have for every p G [1, 00] 

\\e(t)\\ LP <\\e°\\ LP + f \\f(r)\\ LP dr. 
Jo 

We intend to prove the following smoothing effect. 

Theorem 4.2. Let v be a smooth divergence- free vector field o/M d with vorticity u>. 
Then, for every p £ [1, 00 [ there exists a constant C such that 

sup 2«\\A q e\\ L i LP <C\\6 \\ LP + C\\e°\\ L °o\\uj\\ LlLP , 
for every smooth solution 6 of (TD) with f = 0. 



Proof. We start with localizing in frequencies the equation: for q > — 1 we set 
dtd q + v-V6 q + \B\0 q = -{A q ,v V]0. 



9 q := A q 6. Then 
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Recall that 9 q is real function since the functions involved in the dyadic partition of 
the unity are radial. Then multiplying the above equation by \6 q \ p ~ 2 9 q , integrating 
by parts and using Holder inequalities we get 

/ mo g )\o q \ p - 2 e q dx < IIMl^IKa^-vmIlp. 

P dt Jr 2 
Recall from [9] the following generalized Bernstein inequality 

c2 q \\e q \\ p LP < [ (\B\e q )\e q \p-%dx, 

Jr 2 

where c depends on p. Inserting this estimate in the previous one we obtain 

-j^Wlp + ^\\e q \\ p LP < we^WiA^v v]e\\ LP . 

Thus we find 

j t \\e q \\ LP + c2^\\e q \\ LP <\\[A q ,v-v}e\\ LP . 

To estimate the right hand-side, we shall use the following lemma (see the appendix 
for the proof of this lemma) . 

Lemma 4.3. Let v be a smooth divergence-free vector field and 8 be a smooth scalar 
function. Then, for all p £ [l,oo] and q > —1, 

\\[A q ,v .V]6\\ LP <\\Vv\\ LP \\8\\ BOca oo . 
Combined with (|4.ip this lemma yields 

j t (e ct2q \\e q (t)\\ LP ) < e^wvvm^wem^ 

< e ct2q \\uj(t)\\ LP \\0 o \\ Lao . 

To get the last line, we have used the conservation of the L°° norm of 9 and the 
classical fact 

I|Vv||lp < \\u\\lp VpG]l,+oo[. 
Integrating the differential inequality we get 

P q (t)\\ LP < \\0 q \\Lpe- ct2 '' + \\6°\\ L ~ [' e~ c ^ 2q Mr)\\ LP dr. 

Jo 

Integrating in time yields finally 

^ \K\\Lr + \\0°h°° [ MT)\\ LP dT 
JO 

< \\0°\\lp + \\0°\\l°° [ \\cj(r)\\ LP dT, 
Jo 

which is the desired result. □ 

Let us now move to the last part of this section which deals with some logarithmic 
estimates generalizing the results of \2Q\ [16] . First we recall the following result of 
propagation of Besov regularities. 
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Proposition 4.4. Let (p, r) G [1, oo] , s &] — 1, 1[ and 9 a smooth solution of (TD). 
Then we have 

||% rB% <e CV ^(\\e°\\ Blr+ l\~ CV ^\\f(r)\\ B ^dr), 

' •/ 

where V(t) = \\Vv\\ L i Lao . 

The proof of this result is omitted here and it can be done similarly to the inviscid 
case [8] , using especially Proposition 14.11 

Now we will show that for the index regularity s = we can obtain a better estimate 
with a linear growth on Lipschitz norm of the velocity. 

Theorem 4.5. There exists C > such that if k > 0, p € [1, oo] and 6 a solution 
of 

(d t + v-V + K\D\)0 = f, 

then we have 

*t 



T oo dO 



< c(\\e°\\ B o A + ||/|| L i s o ;1 ) (i + / ||v«(T)|| L ~dT) • 



Proof. We mention that the result is first proved in [26] for the case k = by 
using the special structure of the transport equation. In [T7] the first two authors 
generalized Vishik's result for a transport-diffusion equation where the dissipation 
term has the form — kA9. The method described in [T7| can be easily adapted here 
for our model. 

Let q £ NU{— 1} and denote by 9 q the unique global solution of the initial value 
problem 



(4.2) 



dtd q + v-V9 q + \B\9 q = Agf, 

6q\t=0 = A Q 6 °- 



Using Proposition 14.41 with s = ±i we get 

fe q \\ h < (|| vi| ±h +\\A qf \\ LlB±h )e cv «\ 

where V(t) = \\Vv\\ L i Loo . Combined with the definition of Besov spaces this yields 
for j,q > -1 

(4.3) \\Aj9 q \\ L ™ L p < 2-al*-«1(||A,0V + \\A q f\\ LlLP )e cv ^. 
By linearity and again the definition of Besov spaces we have 

(4.4) II^IIloob^ < E \\Aj0qh?LP+ E \\ A ^lh^LP, 

\j-q\>N \j~q\<N 

where TV G N is to be chosen later. To deal with the first sum we use (|4.3p 
£ W^Mifi* Z 2 ~ N/2 £ (l|A^ ||^ + ||A g /|| LtliP )e^W 

\j-q\>N q>-l 

< 2-^(||^|| B o i + ||/|U tlB o i ) e ^). 
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We now turn to the second sum in the right-hand side of (14. 4p . 
It is clear that 

\j-q\<N \j~q\<N 

Applying Proposition 14. II to the system (|4.2p yields 
It follows that 

Y, IIAAIUr^^^OKll^ + ll/llLiBoJ. 

\j-q\<N 

The outcome is the following 



Choosing 



. 1 <(rb ; „ + ii/ii i ,.Bo 1 )(2-'"v"»+iv). 

2C7F(t)- 



N 



log 2 



+ 1, 



we get the desired result. 



□ 



5. Proof of Theorem 11.11 
Throughout this section we use the notation to denote any function of the form 

$k(t) = C exp(...exp(C t)...), 

k times 

where Co depends on the involved norms of the initial data and its value may vary 
from line to line up to some absolute constants. We will make an intensive use 
(without mentionning it) of the following trivial facts 

/ $jfc(r)dr < $ fc (t) and exp( I ^ k (r)dr) < * fc+1 (t). 
io Jo 

The proof of Theorem 11.11 will be done in several steps. The first one deals with 

some a priori estimates for the equations (11 . 2[> . In the second one we prove the 

uniqueness part. Finally, we will discuss the construction of the solutions at the end 

of this section. 

5.1. a priori estimates. As we will see the important quantities to bound for all 
times are ||cj(t)||^oo and ||Vv(t)||£oc. It seems that for subcritical regularities like 
for example H s ,s > 2 or more generally Bp r , s > 1 + | we need only to bound 
the quantity ||<9i#|| L i L oo, which in turn controls ||cj(i)||£,oo, due to Brezis-Galloiiet 
logarithmic estimate, see for example [23] , Even though these quantities seem to 
be less regular than HVuH^i^oo, it is not at all clear how to estimate them without 
involving the latter quantity. 

When we deal with critical regularities which is our subject here one needs to bound 
the Lipschitz norm of the velocity and this will require some refinement analysis, 
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especially Theorem 14.51 seems to be very crucial. To obtain a Lipschitz bound we 
will proceed in several steps: one of the main step is to give an L°°-bound of the 
vorticity but due to some technical difficulties related to Riesz transforms this will 
be not done in a straight way. We prove before an LP estimate with 2 < p < oo and 
this allows us to bound the vorticity in L°° . 

We start with recalling the L p estimate for the temperature function. It is a direct 
consequence of Proposition 14. 1[ 



\m)\\ L p < 11 



Proposition 5.1. Let (v,9) be a smooth solution of (|1.2j) . then for all p G [l,+oo] 

\lp- 



We intend now to bound the LP-norm of the vorticity and to describe a smoothing 
effect for the temperature. 

Proposition 5.2. Lf uj° G LP and 6° G LP n L°° with p g]2, oo[, then 

IKt)||xp + ||%i B i <*i(*)- 

t p,oo 

Proof. Applying Riesz transform TZ to the temperature equation we get 

(5.1) d t K6 + v ■ VTZ6 + \D\ne = -[TZ,v V}9. 
Since \D\1Z = d\ then the function T := lo + 7Z6 satisfies 

(5.2) d t T + v Vr = -[H,v V]9. 
According to the first part of Theorem 13.31 applied with r = 2 we have 

\\[n,v-v]e\\ B0 <\\vv\\ LP (\\e\\ B0 

Using the classical embedding B® 2 w LP which is true only for p £ [2, oo) 

\\[n,vV]e\\ LP < WVvUpiwew^ + weUp). 

Since div v = then we get from the transport equation ()5.2[) 

||r(t)||£P < ||r°|| LP + f \\[n,vV}9(r)\\ LP dr. 
Jo 

Putting together the last two estimates we get 

m)\\LP < lir ^ + 1*1^)11^(11^)1^ + 11^)^ 

< ||a; || LP + \\e°\\ LP + jf* 11^(^11^(11^(^11^ + \\6°\\ LP )dT. 

We have used here the Calderon-Zygmund estimates: for p G (1, oo) 

[|Vu||iP < C||o;|| L p and ||^°|| LP < C\\0°\\ LP . 
On the other hand, from the continuity of the Riesz transform and Proposition 15.11 

Mt)\\ L P < ||r(t)|| LP + ||^|| L p 
< ||r(t)|| LP + ||0°|| LP . 
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This leads to 

\Ht)hr < \\u°\\lp + \\9°\\ LP + f IKt)|| L p(||0(t)|| b o + \\9°\\ LP )dr. 

JO ' 

According to Gronwall lemma we get 

(5.3) \Ht)\\ LP <C e c ^e CmL ' B °^. 
Let TV € N, then we have by Bernstein inequalities and Proposition 

PWlibo 2 < \\s N e\\ LlB0 +\\(id-s N )e\\ L}B0 

t 00,2 t 00,2. t 00,1 

t\\9 Q \\ L ~y/N + Y,\\^A\l}l™ 



< 



q>N 



< v / iV||^ ||^t+^2^||A^|| L i LP . 

q>N 

Using Theorem 14.21 and p > 2 we obtain 

£ 2«l\\A q 9\\ LlLP < Yl ^^(KII^ + KIU- fMr)\\ LP dr 

>N-1 q>N-l J ° 



q>N- 



Thus, we get 



< II0V + 2 



M-i+l)ii fl o 



||0 U || L oo / \\u{r)\\ LP dT. 




LiBL <V^||0°||W + ||0V + 2 



JV(-l+|),mO 



We choose N as follows 



'l|0 u |U~ / Mr)\\ L pdT. 




N 



log(e + / \\uj{r)\\ LP dT 



+ 1. 



(l-2/p)log2 
Then it follows 

H^llijB^ < KHz-rU* + ||e°||L«-tlog3 (e+ f Mr)\\ LP dr). 

J 

Combining this estimate with (|5.3p we get 

II^IIl^o 2 < KlU-ni* + \\e \\ L ~>t \og l 2 (e + Co e Cot e CmL t B ^ 

t oo,2 y 

< coa+^+cii^ii^tii^ . 

Thus we get for every t £ R+ 

II^IIlk, <c (i + t 2 ). 

It follows from (|5.3p 

(5.4) II^WIIlp < ®i(t). 
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Applying Theorem 14.21 and (15. 4h we get 

(5.5) 2«\\A q 9\\ L i LP < *!(t), Vq G N 

and thus 

||e|| £ i Bl < *i(t). 



This ends the proof of Proposition 15.21 

Remark 5.3. It is not hard to see that from (|5.5|) one can obtain that for every s < 1 
(5-6) IMIijBf, < II^IIZ}Bi < $ iW- 

£ Pi -L L -'t p,oo 

Combined with Bernstein inequalities and the fact that p > 2 this yields 
(5-7) MlIb^ < 

for every e < 1 — |. 

□ 

We aim now at giving an L°°-bound of the vorticity. 

Proposition 5.4. Let (v,9) be a smooth solution of (jl.2p suc/i thatuj ,8 E L p r\L°° 
and 1Z6° £ L°° , with 2 < p < oo. Then we have 

(5.8) Mt)h°° + \\ne(t)\\ L °o <* 2 (t) 

and 

(5.9) < ^s(t). 
Proo/. 

• Proof of (|5.8p . By using the maximum principle for the transport equation (|5.2p . 
we get 

||r(t)||L« < ||r°|| L oc + f \\[k,v v]d(T)\\ L °odT. 

Jo 

Since the function 1Z0 satisfies the equation 

(5.10) (dt + v-V + \D\)Ke = -[K,v-V]8, 
we get by using Proposition 14.11 for p = oo that 

\\Ke(t)\\ L - < \\K0(t)\\ LOO + I \\[n,vV)e{r)\\ L oodT. 



o 



Combining the last two estimates yields 

||r(t)IU~ + \\ne{t)\\ LOO < ||r°|| L oo + ||^ || LOO + 2 / \\\ii,vV\6{T)\\ L ~>dn- 

Jo 

< C + [ \\[n,vV]e(r)\\ B0 dr. 
Jo ° c ' 1 
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It follows from the second estimate of Theorem 13.31 and Proposition 15.21 

\\u(t)\\ La o + \\Ke(t)\\ L «, < C + /V(T)IU-nL P (||0(r)|| B;ol + ||0(r)|| iP )dT 

Jo 

< C + \\w\\L°°Lp(\\8\\L}B e + t \\^°\\Lp) 

L ^ t OO, 1 ' 

+ ril^(T)||L-(||e(T)b Ll + ||0 O || LP )dr. 
Jo 

Let < e < 1 — | then using (|5.7p we get 

IK*)|U« + WM(t)h» < m) + f Mt)\\ l ~ {\\e(r)\\ Blai + \\e°\\ LP )dr. 



o 

Therefore we obtain by the Gronwall lemma and a new use of (|5.7p that 

||w(t)||L<» + \\HB(t)\\ L «, < $ 2 (t). 

• Proof of (ESI. 



Let iV € N to be chosen later. Using the fact that ||A 9 v||£°° ~ 2 9 || A 9 u;||l°o , we 
then have 

IK*)||i~ < \\ X (2- N \B\)v(t)\\ L ^ + 2-"||A^(t)|| L o 

q>~N 

< \\ X (2- N \B\)v(t)\\ L ~ +2 N \\u;(t)\\ L ~. 
Applying the frequency localizing operator to the velocity equation we get 

X {2- N \B\)v = x(2- N \B\)v + [ t Vx(2- N \B\)e(T)dr+ f ' V X {2- N \B\)dw(v®v){T)dT. 

Jo Jo 

where V stands for Leray projector. From Bernstein inequalities, C alder on- Zygmund 

estimate and the uniform boundness of x(2 _Ar |D|) we get 

/ < ||x(2- JV |D|)^(r)|| i0 o ( ir < 2~ N l f\\e{r)\\ LP dT 
Jo Jo 

< t\\o°\\ LP . 

Using Proposition I3.U (2) we find 

f \\Vx{2- N \D\)d\v{v ® v){T)\\ L -dT < 2- N f \\v{r)t Laa dr. 
Jo Jo 

The outcome is 

< \\v°\\l°° + t\\e \\ L v + 2~* f \Hr)f Loo dT + 2^11^ 

Jo 

f \\v(T)\\l^dT + 2 N <5> 2 (t) 
Jo 



< 2 



ft 

-N 
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Choosing judiciously N we find 

\\v(t)\\L°° <* 2 (*)(l + (jf Hr^Udr) 
From Gronwall lemma we get 

\\v{t)\\ L °° < $3$). 



□ 



Now we will describe the last part of the a priori estimates. Following the program 
exposed in the beginning, the aim is to get estimates on Vi>. 

Proposition 5.5. Let (v,8) be a smooth solution of (II. 2j) and u°,9° £ B^ ol D LP 
with p e]2, oo[. Then 

\\0(t)\\ B o +|Kt)|| B o +\\v(t)\\Bl , <*3(*). 



Proof. By using the logarithmic estimates of Theorem 14.51 for the equations (|5.1 
and (|5.10p . we obtain 



(5.11) ||r(t)|| B o +||^(t)|| B o <(c + [n,v-v}9 LlB0 )(i + \\Vv\\ L]L , 

oo,l oo.l \ 11 IIL i D co,l/ V * 

Thanks to Theorem 13.31 Propositions 15.41 15.21 and (|5.7p we get 

l|[^^V]0|Li B o < /(IKr)|U» + IKr)M(||(?(r)||B^ 1 + ||fl(r)MdT 

< *2(t). 

By easy computations we get 

||Vv||l°° < ||VA_n;||x,oo + || AgV-^ll^co 

<?GN 

^ [Mix* + ^ ||A g a;||z,°° 

(5.12) < + Mt)|| B o 

oo , 1 

Putting together (|5.1ip and (|5.12p leads to 

\Ht)\\ B o < \\T(t)\\ B o + \\K6(t)\\ B o <* 2 (t)(l+ / HrJIIso dr). 

OO 5 1 OO 5 J. OO 5 X \^ q OO ? 1 J 

Thus we obtain from Gronwall inequality 

(5.13) |Ki)|| B o <$ 3 (i)- 

oo,l 

Coming back to (|5.12p we get 

||V«(t)||ioo < $ 3 (t). 

Let us move to the estimate of t> in the space B^ 1 . By definition we have 

\Ht)\\ B i <\\v(t)\\ L <~ + \\u(t)\\ B0 

oo.l oo,l 
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Combined with (|5.9f) and (j5. 13[) this yields 

\\v(t)\\Bi , < *3(t) 

00,1 

as claimed. To estimate ||0(£)||bo we use Theorem 14.51 and the Lipschitz estimate 
of the velocity. 

11^)11^ £ H^Hb^X 1 + / t ||Vu(r)||L»dr) 
< * 3 (*). 

The proof of Proposition 15.51 is now achieved. □ 

Remark 5.6. Notice that the a priori estimates above imply that V# G L°°). 
Indeed, we can establish the following estimate by combining the smoothing effect 
of the temperature equation with the logarithmic estimate described by Proposition 

ii^ii^^cii^ii^a + iiv.iii^). 

Although we expect to use this estimate for the uniqueness part, it seems that this 
is not sufficient for our purpose and some technical problems arise due to the fact 
that Riesz transforms do not map continuously L°° to itself. The crucial information 
that we need for the uniqueness is 9 S L 1 1 oc (M+; -Bp j00 ). 

5.2. Uniqueness. We now show that the Boussinesq system (jl.2j) has a unique 
solution in the class 

x T = {L^B^ n 4i4,i) x ( l t L p n 4^,oo), 2 < P < 00. 

Let (v 1 , 1 ) and (v 2 , 9 2 ) two solutions of ()1.2j) belonging to the space Xt, and denote 

v = v 2 -v 1 , = e 2_ e i^ 

Then we have the equations 

dtv + v 2 • X7v = —X7p — v • Vv 1 + 6e2 
< d t 9 + v 2 -V9 + \D\9 = -v -V9 1 
v {t=0 = v°, 9 {t=0 = 8°. 

According to Theorem 14.51 we have 

Mt)\\ B o , <^ + V 1 (t))(\\v°\\ B o +\\Vp\\ L i B o +\\v-Vv 1 \\ L i B o +\\9\\ L i B o ), 

oo,l \ oo,l t oo,l t oo,l t oo,l / 

with V\{t) = HVv 1 ]!^!^. A straightforward calculus using the incompressibility of 
the flows gives 

Vp = -VA _1 div (v • V(V +v 2 )) + VA" 1 ^ 
= I + 11. 

To estimate the first term of the RHS we use the definition 

||I|| B o < IKVA-Miv )div A_i(u ® (v 1 + v 2 ))\\ L oo + \\ v ■ V(v l + v 2 )\\ Bl 

oo ,1 00 , 1 
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From Proposition 13.11 (2) and Besov embeddings we have 

||(VA- 1 div)div A-!(v ® (v 1 + v 2 ))\\ L oo < \\v ® (v 1 + v 2 )\\ L oo 

i$ \\v\\b° \\v 1 + v 2 \\ b o 

00,1 00,1 

Using the incompressibility of v and using Bony's decomposition one can easily 
obtain 

Putting together these estimates yields 

(5-14) ||I|| B o < \\v\\ B o Wv 1 +v 2 \\ B i 

00,1 00,1 00,1 

Let us now show how to estimate the second term II. By using Besov embeddings 
and C alder on- Zygmund estimate we get 

||II|| B o < IIVA- 1 ^ 



00 , 1 

< 



Combining this estimate with (|5.14p yields 

\Ht)\\ B ^ <z cv{t) {\\A\ B ^ + [ 11^)11^ [l+ll^yjWH^Jdr) 

(5.15) + e cv ^\\9\\ 2 , 

where V(t) := IK* 1 ,* 2 )^^. 

We have now to estimate lid II 2 . By applying A to the equation of 9 and arguing 
similarly to the proof of Theorem 14.21 we obtain 



\e° g \\ LP + I e- c2q ^\\A q {vVe l ){r)\\ L pdr 
Jo 

+ f e^ 2q ^\\[v 2 -V^ q ]e{r)\\ LP dr. 
Jo 



Remark, first, that an obvious Holder inequality yields that for every e £ [0, 1] there 
exists an absolute constant C such that 

f e- cr2q dT < Ct £ 2~ qi - l - £ \ Vt > 0. 
Jo 

Using this fact and integrating in time 

9 9 fllfl II , < 9 9 ( _1+ f)|lfl°ll r 
1 p \\tiq\\L]Lv ~ 1 v \Wq\\LP 

+ t e 2 q{ ~ l+e+ l ) J (||A g ( V -V0 1 )(r)|| LP + ||[ ? ; 2 -V,A g ]e(T)|| L p)dr 
(5.16) = 2 9( - 1+ l ) ||^||^+I g (t) + II g (t). 
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Using Bony's decomposition we get easily 

HA^.v^Xt)!!^ < \\v(t)\\»- E n^o\t)\\L* 

+ nm\u~ E \\wHt)\\i*- 

Integrating in time we get 

h(t) < t e ||HlL r ^2^ 1+ l +e) (k| + l)||0 1 || ZlBl 

+ t £ \\v\\ LrL ~2 q{1 P +£) W^WlIlp 
j>q-<i 

(5-17) < t e |b||L r L-2 9( - 1+ l +£ )(|g| + l)||0 1 || ZlB1 . 

To estimate the term 11^ we use the following classical commutator (2/p < 1) 
\\[v 2 .V,A q ]9\\ LP <2^\\Vv 2 \\ L ~\\e\\ a . 

Thus we obtain, 

(5.18) lI q (t)<t s 2^- 1+ ^\\Vv 2 \\ LrL o 



T 1 R p 

t PA 

2 
V 

summing ()5. 16j) and using (|5.17p and (|5.18p we get 



We choose e > such that — 1 + - + e < 0, which is possible since p > 2. Then 



II0II 2 £ \\d°\\L P +t £ \\v\\ LrL o,\\9 1 \\ llBl +t £ ||V7; 2 || LrLO o||^|| a . 
For small i E [0, 5] one can obtain 

t p,l 

Plugging this estimate into ([5.150 we find 

IML~b° , <e cyW (lb°|| B o + ||0°|| LP +*HI o +t e ||t;|U r £<»||e 1 || riBl 

fi oo.l \ oo,l t oo,l 1 t p, 

If 5 is sufficiently small then we get for t E [0, 5] 

(5-19) |M| irB o ^KIIbo ^HflV- 

t 00 , 1 oo,l 

This gives in turn 

(5.20) \\9\\ a <|| W || B o i + ||0 o |U,. 

rlnP oo,l 

This gives in particular the uniqueness on [0, 5]. An iteration of this argument yields 
the uniqueness in [0, T]. 



(B„) 
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5.3. Existence. We consider the following system 

d t v n + v n ■ Vv n + Vp n = 9 n e 2 

d t e n + v n -ve n + \B\e n = o 

divv n = 

V n \t=0 = S n V°, 9 n \t=0 = S n 9° 

By using the same method as |16j we can prove that this system has a unique local 
smooth solution (v n ,6 n ). The global existence of these solutions is governed by the 
following criterion: we can push the construction beyond the time T if the quantity 
|| ViVi || £i £oo is finite. Now from the a priori estimates the Lipschitz norm can not 
blow up in finite time and then the solution (v n ,6 n ) is globally defined. Once again 
from the a priori estimates we have for 2 < p < oo 

Then it follows that up to an extraction the sequence (v n , 9 n ) is weakly convergent 
to (v,9) belonging to LfB^ x L5?(B°, }1 n L p ), with uj e LfL v . For (n,m) E N 2 
we set v ntm = v n — v m and Q n ,m = 9 n — @m then according to the estimate (|5.19p 
and (|fT2o| we get for T = 5 

\\ v n,m\\ x \\9n,m 

L T B pl 

2 

This shows that (v n , 9 n ) is a Cauchy sequence in the Banach space L^B^^ x L\B^ X 
and then it converges strongly to (v,9). This allows to pass to the limit in the system 
(B n ) and then we get that (v, 9) is a solution of the Boussinesq system (jl.2p . 



Appendix: proof of Lemma [4751 

We have from Bony's decomposition 

[A q ,vV]9 = V]A j 9 + Y [A q ,AjV ■ V]-%_i0 

|i-?|<4 li-g|<4 

+ Y [\,AjV-V]Aj6 

i><?-4 



I, + II, + III,. 



Observe first that 



l q = Y h q *(Sj-iv ■ VAj9) - Sj-iv ■ (h q * VAj9) 
b-g|<4 

where h q (£) = ip{2~ q £ > ). Thus, Lemma 13.21 and Bernstein inequalities yield 
\\W\lp < Y \\ x K\\LAWSj-iv\\ LP \\V Aj9\\ L °° 

b'-3l<4 

< ||V«|M|/*,|U* Y V~ q \\Aj h°° 

b'-5l<4 

< IIVv||lp||0||bo 
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To estimate the second term we use once again Lemma 13.21 

IIUJlp £ Y 2-"||A,V V || L p||V5 i _ 1 ^|| L - 

b-<7l<4 

< \\Vv\\ LP Y 2 fc ^||A^|| L c 



\j-q\<i 

k<j-2 



< HVul 



Let us now move to the remainder term. We separate it into two terms: high 
frequencies and low frequencies. 

Ill, = Y [A^A^A^ + tA^A-xu. V]A_x0 
:= IIlJ+III^. 

For the first term we don't need to use the structure of the commutator. We estimate 
separately each term of the commutator by using Bernstein inequalities. 

Iimjllw < Y 29 \\ A j v hA\^h^ 

j>q-4 

< I|Vu||xp Y 2"" J ||A^IU- 

J>9-4 

< UVullipllflllflo,^. 

For the second term we use Lemma 13.21 combined with Bernstein inequalities. 

PUgllLp < ||VA_ii;||ip||VA_i0||£cc 
< ||Vt;| 
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